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THE  CALCULA.TION  OF  A  THEKMAL  BOUNDARY  LAYER 
IN  THE  FLCW  OF  A  COMPRESSIBLE  GAS 


By 

L.  M.  Zyslna-Molozhen 

A  seml-emplrlcal  approximation  method  which  permits  calciilatlon 
with  adequate  accuracy  of  the  laminar^  transient  and  turbulent  regions 
of  a  thermal  boundary  layer  during  flow  of  compressible  gas  over  its 
surface,  is  examined  in  this  aurtlcle. 

The  solution  of  many  technical  problems  is  often  associated  with 
the  necessity  of  calculating  the  heat  exchange  of  a  surface  with  the 
compressed  gas  flowing  over  it.  The  boundary  layer  appearing  imder  such 
circumstances  and  dependent  on  the  flow  conditions  and  the  nature  of 
velocity  distribution  along  the  surface,  may  either  be  laminar,  transient 
or  turbiilent  on  the  overall  surface;  or  it  may  be  laminar,  then  transient, 
and  then  turbulent  on  perts  of  the  surface. 

In  this  article  a  seml-empirlcsil  approximation  method  is  proposed, 
which  permits  calculation  with  sufficient  accuracy  of  all  three  thermal 
regions  of  a  boundary  layer  appearing  during  the  flow  of  compressible  gas 
over  its  surface.  The  essence  of  the  method  Is  as  follows. 

Let  us  examine  a  plane  flow  of  a  compressible  gas.  It  is  known  that, 
in  this  case,  the  integral  relationship  of  energy  in  Dorodnitsln  variables 
appears  as: 
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(sign  «  pertains  to  the  parameters  of  drag) . 


It  may  be  vritten  as 


d  ±.  d 

-jfO- •?)-' -5T  + 


Dit  0 

dTln  ’ 


(5) 


where 
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Using  relationship  (5),  equation  (1)  may  be  reduced  to  appear  as: 


(1  -«2) 


Nu 

PrRej 


Introducing  the  para¬ 
meters: 


Nu* 


(8) 

(9) 
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let  US  assume  that -Qiey,  changing  along  the  streamlined  surface,  uniquely 

determine  all  characteristics  of  a  thermal  boundary  layer*  Let  us 

also  assume  that  In  expression  (8)  the  effect  of  a  longitudinal  gradient 

U*  8**  ^ 

of  pressure  Is  characterized  by  conplex  and  Q-  Is  not  dependent 

on  ~  and  only  determines  the  effect  of  Reynolds  number*  In  this  case, 
value  vUl  be  the  same  for  the  flow  around  the  profile  as  well  as  along 
the  lamina*  We  determine  function  dip  by  data  on  the  heat  exchange  of  the 
lamina* 

By  examining  the  test  data  (Fig  1)  It  Is  apparent  that  dtirlng  deviation 
of  the  physical  constants  toward  the  parameters  of  drag*,  the  formulas  for 
calculating  heat  exchange  retain  the  same  form  as  in  the  case  of  an 
incompressible  flow* 

For  the  transient  region  of  a  boundary  layer  the  influence  of 
compressibility  reflects  on  the  coordinates  of  the  start  and  the  finish  of 

transition  (■»■■«.  and  the  development  of  the  transition  process,  after  its 
appearance,  occvu^  in  such  a  manner  that  the  linos  Nu,  =-  A'(Re,)  remain 

parallel  to  each  other  in  the  whole  investigated  diapason  of  the  change  in  the 
number  M* 

Evidently  these  curves  can  be  approximated  by  series: 

Nu,  =  BRe;,  (10) 

where  values  B  and  n  are  different,  but  are  invariable  for  each  regime  of  flow 
in  the  boundary  layer*  For  a  transient  region  coefficient  B,  as  the  teste  show, 

^During  the  calculations  all  physical  constants  deviate  toward  the  tenpera- 
ture  of  drag. 
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is  a  variable  value,  changing  uith  the  change  din  the  value  of  Rejn^olda  nuniber, 
vhich  corresponds  with  the  initial  point  of  transition.  This  value  Is 
deterndned  by  the  initial  turbulence  of  tho  flow,  by  the  value  of  the 
temperature  factor  and  the  number  M«,  i«e.,  value  B  retains  a  constant  value 
within  the  limits  of  each  concrete  test,  but  can  change  with  the  chutge  of  the 
Initial  and  operational  conditions  of  the  process. 

Substituting  expression  (10)  into  the  integral  energy  relationship  for 
a  lamina: 


T* 


Nu, 


(11) 


and  using  formula  (9),  we  can  obtain  the  following  expression  for  Qipt 


[Pr 
- 1 


'  T'  *  I 


k  ijR-fl  RCf  ’ 


.  _ 1  — a 

where  ^  *•  * 


We  determine  function  X  in  such  a  manner  that  the  following  may  be 
acceptable  for  the  lamina: 


r  T:  J.1-<m+i) 


(13) 


Then  the  formula  for  Qj,  Is  ob¬ 
tained: 


Cr.»^Re”  , 


(W) 


which  is  conpletely  analogous  to  the  f omuls  for  an  Ineompresslble  flow  /!/. 

k 


{ 
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Fig.l,  Dependence  of  “  N  (Rejj)  at  different 
values  of  the  nuroer  M:  2-0«l4i»} 

3-0.52J  lt-0.73j  5-1.07J  6-1.1(35  I-Nu^  -0.297 
Re^  1/2-  laminar  boundary  layer;  H-  Nu*  ■  0,0233 
Rex  U/3-  turbulent;  VT^*  ~  1, 

By  introducing  relationships  (8),  (9),  and  (lU)  into  equation  (7)  and 
performing  siiqple  calculations,  it  is  possible  to  reduce  it  tb  the  following 
form: 

■  (.6) 

For  incoiqpressible  flow  it  was  experimentally  shofwn  /2/  that  the  function, 
placed  in  brackets 

fr-(m+ l)x-2/r  (16) 
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may,  with  sufficient  accuracy,  be  expressed  by  relationship: 

fr-a-2/r.  (17) 

With  this  it  was  found  that  value  c:  is  not  dependent  on  the  gradient 
of  pressure,  but  is  determined  only  by  the  regime  of  flow  in  the  boundary 
layer. 

let  us  assume  that  an  analogous  relationship  nay  also  be  written  for  a 
compressed -flow.  If  a  is  not  dependent  on  the  longitudinal  gradient  of  the 
pressure,  then  it  evidently  should  have  the  same  value  for  the  flow  along 
the  profile,  as  well  as  for  the  flow  along  the  lamina.  Then,  using  expression 
(13),  OL  nay  be  determined  by  formula 


0“  (/n+  1) 


(18) 


where  is  value  cc^  for  the  undisturbed  flow, 

VJhen  examining  the  flow,  characterized  by  conditions  T ®  “  const, 
it  is  evident  that  for  each  regime  of  the  flow  a  will  be  a  constant  value 
and  dependent  only  on  the  regime  of  the  flow. 

In  this  case  equation  (l^)  is  easily  integrated  and  allows  to  determine 
6*^  for  the  thermal  boundary  layer: 


I  m  .V  ^ 
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Here  corresponds  to  value  of  coefficient  A,  in  formula  (lU)  for 
the  laminar  boundary  layer  if  the  calculation  is  carried  out  for  a  transient 
region,  and  corresponds  to  value  A  for  a  transient  region  if  the  calculation 
is  carried  out  for  a  turbulent  section  of  the  boundary  layer.  Accordingly, 
Re’*  »  ^0  ^Tn  is  determined  for  coordinates  of  the  start  of  transition  from 

‘"'T'h  — r— 

^0 

calculation  of  the  laminar  section  in  the  first  ease,  and  for  coordinates 
of  the  finish  of  transition  from  calculation  of  the  transient  region  in  the 
second  case.  For  convenience  and  speed  up  of  the  calculation,  values  ct  for  ’ 
each  regime  of  flow  in  the  bo^ondary  layer  may  be  calculated  by  formula  (18) 
in  the  required  diapason  of  changes  of  paramteirs  and  be  presented  in  the 
form  of  curves  /Fig  2  /, 

Using  formulas  (19),  (9)»  (l6)  -  (18)  an  e:qpressicn  may  be  obtained  for 
calculation  of  local  values  of  the  heat  exchange  coefficient 


By  conparing  (20)  with  a  corresponding  formula  for  an  tncoi)5)ressed 
flow(U)  it  is  discovered  that, during  deviation  of  all  physical  constants 
toward  the  ten^erature  of  drag  of  flow  ,  a  formal  analogy  exists  between 
the  appearance  of  the  formula  for  calculation  of  the  heat  exchange  intensity 

I 

in  a  con^ressed  flow  Nu^  ,  and  an  incoit^ressed  flow  Nu^  •  By  cosparing  these 
formulas  a  relationship  is  easily  obtained 
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Fig  2,  •  Dependence  of  a  ■  O'(X)  for 
a  laminar  region  of  the  boundary 
layer  (a);  for  a  transient  region  (b) 
and  for  a  turb\ilent  region  (c)  at 
different  T yT*  ; 

/-0,5;  2—0,0  0,7;  <1^0.8;  5-0,0;  5—1.0;  7— 

I,?;  5-K6;  /O— 1,8;  //— 2.0 


It  is  seen  from  this  formula  that 
in  the  instance  of  a  flow  over  the 
lamina  with  velocity  on  the  outer 
boundary  layer  equaling  the  velocity 
of  the  Inflowing  stream  and  corresponds 


the  above  described  representation  of  the  formulas,  the  following  relationship 
is  valid: 

0 

Nu,  «=  NujI-  (22) 

As  was  already  mentioned,  this  formal  relationship  corresponds  well  with  the 
test  data. 

During  calculation  the  coordinate  points  of  the  start  of  transition 
may  be  determined  by  Dorodnltsin  -  Loitsinskii  method  /3/«  For  an 
approximate  determination  of  the  transient  region,  the  folloiTing  siitple 
consideration  may  be  recommended.  Processing  of  the  tests  shown  In  Fig,  1 
and  analogous  tests  of  other  authors  indicated  that  pararater  r^,  characterizing 
the  relationship  between  the  coordinates  of  the  finish  and  the  start  of 
transition  do  not  depend  on  the  number  H,  This  characteristic  of  parameter 
r^  permits  determination  of  its  value  along  eitpirical  curves  obtained  for 
an  inconpressed  flow  and  presented  in  article  /V* 

DESIGNATIONS 

0^  is  the  density  corresponding  to  parameters  of  drag;  U^,  T^  is  the  velocity 
and  the  temperature  outside  the  boundary  layer;  T^  is  the  tenperature  of  the 
wall. 
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